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Abstract 

We find conditions required to achieve complete population trans- 
fer, via coherent population trapping, from an initial state to a des- 
ignated final state at a designated time in a degenerate n-state atom, 
where transitions are caused by an external interaction. In systems 
with degenerate states there is no time ordering. Analytic expressions 
have been found for transition probabilities in a degenerate n-state 
atom interacting with a strong external field that gives a common 
time dependence to all of the transition matrix elements. Except for 
solving a simple n^^ order equation to determine eigenvalues of dressed 
states, the method is entirely analytic. These expressions may be used 
to control electron populations in degenerate n-state atoms. Examples 
are given for n = 2 and n = 3. 
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1 Introduction 



Analytic descriptions of physical systems are convenient for understanding 
quickly and easily how a system works under various conditions. An example 
is population control in quantum systems, namely transfer of electrons from 
an ensemble of atoms all in the same initial state to specified final states 
within the ensemble. This is used in problems ranging from quantum infor- 
mation m 12] to chemical dynamics E] • Such problems are often studied in 
terms of a single electron in an atom with n discrete states interacting with a 
strong external field 13 El 13 UHl HI]- In this paper analytic expressions 
are presented for the transition probabilities as a function of time in a de- 
generate n-state atom. Once analytic solutions are obtained, corrections due 
to finite energy splittings can be introduced, e.g. in numerical calculations. 
Then real systems with finite differences in the energies of the n states can 
be studied. 

Our degenerate energy approximation is somewhat similar to the rotating 
wave approximation (RWA) E] that has been widely applied to both 2 and 
3-state atomic models. In RWA, however, degenerate atomic states are not 
used. Instead one tunes the frequency of the external field V{t) = Vo cos(a't) 
to the frequency difference of two non-degenerate states so that the detuning 
parameter A = h{ijj — UJ12) tends to zero, where UJ12 = {Ei — E-2)/Ti. Thus in 
RWA an initial state of an atom plus one photon is degenerate in energy with 
the final state of the atom. An advantage of using degenerate atomic states 
is that one is not restricted to external interactions with frequencies close to 
the transition frequency. Thus uj can be used for control, e.g. to vary the 
duration of time that the transferred population remains in the designated 
state, or to reduce the population leakage that occurs when the energy states 
are not fully degenerate. On the other hand RWA has the advantage that 
analytic solutions have been found when the detuning parameter is small, 
but finite. 

2 Theory 

Let us consider an n-state atom interacting with an external field, Vextif^t). 
The total Hamiltonian for this system is H = Hq + Vexti'f^it)- The n eigen- 
states, and corresponding eigen-energies, E^, of Hq are assumed to be 
known. The total wavefunction may be expanded in terms of the known 
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eigenstates, namely, = ai(t)0i + a2(t)02 + ■ ■ ■ + an{t)4>n- With atomic 
units, using i4f = {Hq + V;^t(f, t))^, with ifo0fc = Ek4>k and J(f)*4>kdr = Sjk, 
one then obtains 0, 

n 

iaj{t) = Ejaj{t) + J2Vjk{t)ak{t) , (1) 

A;=l 

where Vjk(t) = J (f)*Vext{r,t)(f)kdr. These equations are exact for an n-state 
atom. 

We assume that the system is degenerate, namely that all the energies, Ej, 
are the same. Since the zero point of energy is arbitrary, one may generally 
set Ej = 0. These conditions give the coupled equations for our degenerate 
n-state system, namely, 

n 

idj{t) = Y,Vjkit)ak{t) . (2) 

k=l 

We additionally require that all of the Vjkit) have the same time dependence, 
e.g. Vjkit) = VjkEofit), where f(t) is a real, but otherwise arbitrary function 
of time. For this work we consider Vjk = Vkj to be real. We use the initial 
conditions ai(0) = 1, and 0^(0) = for j ^ 1. 



2.1 2-state atom 

We first illustrate our scheme for population control in a two state atom. For 
n = 2 Eq(j21) becomes, 

idiit) = \/ii(t)ai(t) + \/i2(t)a2(t) 

la^it) = V2i{t)ai{t) + V22{t)a2{t) . (3) 

Set V^2{t) = V2iit) = V{t), and V.jit) = e,V(t). 

Now seek a solution c(t) = xiai(t) + X2a2(t) such that ic(t) = zV(t) c(t), 
where z is to be determined. Since ai(0) = 1 and 02(0) = 0, we may set 
Xi = 1. Using Eq(jni), one has, 

ic{t) = idiit) + X2 id2{t) 

= eiV{t)ai{t) + V{t)a2{t) + X2 V{t)ai{t) + X2 e2V{t)a2{t) 

= (x2 + ei)\/(t)K(t) + i^^a2(t)] 

X2 + ei 

= zV{t)c{t) . (4) 
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This holds if and only if z = X2 + ei and ^^^^ = X2. This yields two 
roots for X2 and z, namely, X(i^2)2 = ± y^l + (^^^)^ and 21^2 = ± 
^JY+IJ^^^. Using X12 = X(i)2, X22 = X(2)2, xii = X21 = 1, and defining the 
action, A{t) = /q V{t')dt', the corresponding eigenfunctions are, 

ci(t) = e-*^i^W = xnai(t) + Xi2a2(t) 

C2(t) = e-'^^^W =X2iai(t)+X22a2(t) • (5) 

The relation between the Ci(t) and aj(t) may be expressed as Cj(t) = J2] -MijO^jit), 
where A4ij = Xij. This matrix may be inverted to give the population am- 
phtudes, 

ai{t) = ^[x22Cl(t) - Xi2C2(t)] 

a2(t) = ^[-a;2ici(t) + a;iiC2(t)] , (6) 



where A = 2y 1 + (^^-g^)^ is the determinant of Aiij. 
A simple solution occurs when ei = €2 = e, namely, 

a,{t) = e-'^^W ^[e-'^(*) +e+'^W] = e-'^^Wcos(A(t)) 

a2(t) = e-'^^(*) i[e-'^(*) -e+*^(*)] = -ze-'^^(*)sin(A(t)) . (7) 

One may now determine the conditions under which the population of 
state 2 takes on any desired value P2 = Pv = v"^, where < < 1. Note 
that Pi = I — P2. For EqQ this condition is satisfied if a2(t) = sin(y4(t)) = v 
(within an overall phase), i.e., A{t) = sin~^(t;). This can always be satisfied. 
Note that for a given V(t) one may choose the time such that A{t) = sin~^(f ). 
Alternatively to move population to this arbitrary value at a particular time 
to, one may adjust V{t) so that A{to) = sin~^(t>). If ei 7^ 62, it is easily shown 
that P2it) <!/(! + (e2 — ei)^/4) < 1. Thus the diagonal matrix elements of 
Vij{t) prevent complete transfer to the initially unoccupied state when they 
are unequal, and simply contribute an overall phase when they are equal. 

An interesting solution occurs when the ratio x/^ = ^/2, see ^Sl- Then 
Eqs ((Tj) take the form 

i-TT 

ai = cos[— sm{ujt)\ 

71 

a2 = z sin[— sin(u;t)] . (8) 
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Time (in units of the field period, T=27:/a)) 

Figure 1: Analytical solution for x/cj = 7r/2 for a system with two degenerate 
energy states. 

The corresponding occupation probabilities Pi = |aip and P2 = |a2p are 

Pi = cos^[— sin(co't)] 

P2 = sin2[^sin(cut)] . (9) 

For this special value of x/^ the probability is completely transfered back 
and forth between the two states with the period vr/u; as shown in figure ^ 
It can be easily shown that when P2{t) reaches its maximum value at t = 
7i/2uj, a small deviation e from the value 7i/2 in the ratio x/^ reduces P2 
from 1 to 1 — e^. This can be used to estimate how sensitive this approach 
is to fluctuations in the strength of the external field, x- The period of 
oscillation of the occupation probabilities is two times smaller than the period 
of oscillation of the radiation field since Pi = |ajp. 

Complete population transfer occurs in general wherever the action, A{t) = 
lo ^i(t) dr, is an integer multiple of h/A. Each cycle has a "fiat" part where 
the probability value remains very close to the extreme values of or 1 for an 
extended period of time (cf. figure Q). As we can see, the external radiation 
field with the ratio x/^ = can be used to cause periodical population 
inversion of the electronic states (a physical realization of a quantum con- 
trollable system). 

Equations (jH} can be expanded in a power series in the vicinity of the 
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Figure 2: Population leakage in the vicinity of to = T/4. Solid line - analytic 
solution dash-dotted line - polynomial approximation ([TU)). 



point to = T/4, where P2{t) reaches its first maximum. Then P2 can be 
accurately approximated by taking only first few terms in the series, which 
converges rapidly for r < T = 27r/uj. If one chooses the potential in the form 
Vexti'i^yt) = —fE(t), the first, second and third derivatives of the function 
P2 are zero at t = Iq, and 

P2(to + r) = 1 + i^r^ + 0{r') ^ 1 - ^{urf • (10) 

Hence the occupation probability P2 in the vicinity of its maximum can be 
approximated by a Ath degree polynomial. 

The polynomial approximation of Eq. ()10j) tells one how to choose the field 
frequency to provide a desired duration of the populated state. To obtain a 
duration Tg of the populated state with a population leakage less than some 
critical value Per, the field frequency to be used is 




Near r = the error is 0{{ujt)^). This is illustrated in figure El 
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2.2 3-state atom 



Next consider a 3-state atom. The general case for n = 3 can be solved 
by taking V{t) = a~^Vi2{t) = [3-^Vr^{t) = V2^{t) = e-%j{t). Then Eq© 
becomes, 

= eiV{t)ai{t) + aV{t)a2{t) + pV{t)a-:i{t) 
ia^it) = aV{t)ai{t) + e2V{t)a2{t) + V{t)a3{t) 
2a3(t) = pV{t)ai{t) + V{t)a2{t) + esV{t)a:,{t) . (12) 

Try c(t) = Xiai(t) + X2a2{t) + 0:303 (t). Again, Xi = 1, and for convenience 
we set X2 = X and X3 = y. Then, 

ic{t) = iai(t) + X ia2(t) + y ia^it) 

= eiV{t)ai{t) + aV{t)a2{t) + pV{t)a3{t) 
+x[aV{t)ai{t) + e2V{t)a2{t) + V{t)a3{t)] 
+y[(3V{t)a,it) + V{t)a2it) + esVit)as{t)] 

= {ei + ax + py)V{t) 

a + e2X + y p + x + e^y 

x ai(t)H 7^a2(t) H 7^«3(i) • (13) 

^ ei + ax + (3y ^' ei + ax + (3y^'^ ^ ' 

We seek c(t) such that ic{t) = zV{t)c(t). The last line holds if and only if z = 
ei+ax+/3y, {a+e2X+y) / {ei+ax+/3y) = x and (fS+x+esy) / {ei+ax+/3y) = y. 
It can be easily shown that y = (q;(x^ — 1) + (ei — e2)x)/{l — j3x) and 
X = {l3{y'^ ~ 1) + (^1 " ^z)y)/ (1 ~ oiy)- After some algebra this leads to the 
cubic equation, [(a^ — /9^) +a/3(e3 — 62)]^^ + [/3(2 — — /5^) +a(2ei — 62 — 63) + 
/3(ei - €2) (€3 - e2)]x2 + [(2/32 - a2 _ 1) + ct/3(2e2 - ei - €3) + (ei - 62) (ei - e3)]a; + 
[/?(q;2 — 1) — tt(ei — £3)] = 0. This yields three eigenvalues for x (namely Xj) 
and three eigenvalues for y (yj). 

The eigenvalues, {xj} and {yj}, determine eigenvalues for {zj}, and three 
eigenfunctions, Cj{t) = e~*^j^W, where A{t) = jQV{t')dt'. Specifically, cj = 
T,^=iMjiai, where, 

/I xi yi\ 
M=\ 1 X2 y2 \ (14) 
V 1 3:3 ys J 

This matrix, Ai, may be inverted, namely, 

^ X2I/3 - X3y2 xsyi - xiys Xiy2 - X2yi 

y2 -ys ys- yi z/i - z/2 | (15) 

\ X3 — X2 Xi — X3 X2 — Xi 



A 
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where A = det(A^) = Xiy2 + X2y^ + x^yi — Xiy^ — X2yi — x^y2- Now one may 
express the unperturbed state amphtudes, aj(t), in terms of the dressed-state 
amphtudes, Cj{t). Using Cj{t) = e~*^j^(*\ one has ai{t) = Y.^=i M~j^Cj{t) = 
EjLi M^^e-''^^^'\ Specifically, 

ai = ^ ((X22/3 - X32/2)e-^^^^W + (xaz/i - xiys)e-'''^^'^ + {xiy2 - X22/i)e-^^^^«) 

as = ^((x3-X2)e-^^i^W + (xi-X3)e-*^^^W + (x2-xi)e-*^^^W) . (16) 
The transition probabilities are given by the analytic expressions, 

Pi{t) = lai^t)]"^ = ^[{x2y3 - X3y2f + {xsyi - xiysf + {xiy2 - X2yiy 
+2{x2y3 - X3y2)ix3yi - Xiy3) cos((2;i - Z2)A{t)) 
+2(x2?/3 - 3^3^/2) (a;i2/2 - 3:2^1) cos((2;i - Z3)A{t)) 
+2(x3?/i - Xiy3){xiy2 - X2yi) cos{{z2 - Z3)A{t))] , 

P2{t) = \a2{t)? = ^^[{y2-y3f + {yi-yif + {yi-y2f 

+2(^/2 - 1/3) (2/3 - yi) cos{{zi - Z2)A{t)) 
+2(1/2 - 1/3) (1/1 - 1/2) cos((2;i - Z3)A{t)) 

+ 2(?/3 - yi)(yi - y2) C0S{{Z2 - Z3)A{t))] , 
Pi{t) = \a3(t)\^ = -^[{X3 - X2)^ + {Xi - Xsf + {X2 - Xi)'^ 

+2{x3 - X2){xi - X3) cos((zi - Z2)A{t)) 
+2{x3 - X2)(x2 - xi) cos((2;i - Z3)A{t)) 

+2(xi - X3)(X2 - Xi) cos((2;2 - Z3)A{t))] . (17) 

One may now seek the conditions on the external field V{t) such that the 
electron populations Pj(t) = |aj(t)p take desired values. We have been able 
to show ^2] that for n = 3 the electron can be fully transferred to a targeted 
state at an arbitrary time t = to by suitably adjusting the Vij(t). 

2.3 n-state atoms for n > 4 

For a 4-state atom one similarly obtains a 4^^ order equation in x, yielding 
four eigenvalues for Zj and four corresponding eigenf unctions Cj{t) = e~*^j"^(*). 
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For n > 4, analytic solutions to an n order equation do not exist. However, 
one may determine the eigenvalues Xj by the numerical method of successive 
approximations [12]. It is possible to control the occupation probabilities, 
Pkit), for any n by adjusting the magnitudes of the Vij, e.g. a, (3 and V in 
Eq((T2)), and by changing the shape of V{t). 

As will be verified below, there is a general scheme to find analytic so- 
lutions for degenerate n-state atoms ^3]. This scheme is straightforward. 
Seek a solution to Eq© of the form c(t) = xiai{t) + X2a2(t) + ■ — h x„a„(t). 
Since ai(0) = 6ii, one has that xi = 1. Next calculate ic(t) using Eq© 
and require that ic{t) = zV{t) c{t). Here V{t) is a common factor for the 
Vjk{t) terms in Eq©, and 2; is a linear combination of the Xj's, dependent 
on the relative (time independent) strengths of the Vjkit). This leads to an 
^th Qj-^Qj. equation in X2 (or any of the other Xj's {i > 2)), whose roots may 
be denoted by Xj2 (or Xji in general). This yields n eigenvalues, Zj, and n 
eigenfunctions, Cj{t) = e~*^J^W, where A{t) = jQV{t')dt'. This process de- 
termines the matrix elements, J^ij, for Cj(t) = J^J -M-ijajit). Specifically, 
Aiij = Xij. Inverting this relation yields the probability amplitudes for 
the electron population, afc(t) = M^jCj{t) = Y.] M^je-''^^^^\ Using 
cos(a — b) = cos a cos 6 + sin a sin b, one quickly obtains, 

n n 

Pk{t) = \a,it)\' = E E M-,lM~,} cos[iz. - zMit)] • (18) 

i j 

Since the x^/s and Zj^s vary with the Vjk{t), one may seek conditions on 
the matrix elements Vjkit) and on A(to) such that the electron populations 
Pkito) = |afc(^o)P take desired values at t = to. It has been shown [T21II3] that 
complete population transfer occurs in 2-state and 3-state atoms at t = to if 
A(to)/7r = 1/2 in the 2-state atom and A{tQ)/TT = I/V2 in the 3-state atom. 
In addition for the 3-state atom Vi2{t) = and Vis^t) = V23(t). 

We now note that the n-state equations simplify if Vij{t) = 7V23(t) for all 
J > 3, since all 0^(0) = 0. Using V{t) = a-^Vi2 = f^-^Vn = 7"'V^23 = e]Wu, 
Eq(jH) becomes, 

iaiit) = V{t) {ei ai{t) + a a2{t) + {n-2)P asit) ) , (19) 
id2{t) = V{t) {a ai{t) + e2 a2{t) + {n-2)-f asit) ) , 

id^it) = V{t) ( 13 ai{t) + 7 a2{t) + (es + ^^) ag ) . 
Since V{t) is arbitrary at this point, we may set 7 = 1 without loss of 
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generality. Note that these equations differs from the equations for true 
3-state atoms in that Vij{t) ^ Vji(t). 

These equations may be solved using the general scheme described above. 
To find a solution that is mathematically simple, following the 3-state atom 
[T^ we choose /3 = 1 and = e. The value of e may be arbitrarily changed 
by an overall phase transformation of the a^. Taking X2 = x and = y 
to simplify notation, one quickly obtains, x = {a + y)/ {ax + y) and y = 
{{1 + x) + {^^)y) / {ax + y) . This yields cubic equations in x and However, 
it is evident in this case that if a; = — 1 then y = and if x = 1 then y = y± = 



2(— « + — 3± y {a + n — + 2/ {n — 2). Hence there are three eigenvalues 
for X, y and z, namely {xi} = {1,1,-1}, {yi} = {y+,y-,0} and {zi} = 
{y+ + a,y^ + a, —a}. This gives three eigenf unctions, Cj = e~*^j^(*\ where 
A{t) = /o* V{t')dt', which linear combination of the ai{t). Specifically 

/ -y^ y+ {y+ - y-) 



Adij^Cj, where, ^A 



J, ,,..v.xv., .... 2{y+^y^) 



-y- y+ -{y+-y-) | This yields 
V 1 -1 
transition probabilities Pi{t) = |aj(t)p, in accord with Eq (fTK|) . 

From Eq p8|l extrema of Pi{t) occur at t = to when {zi — Z2)A{tQ)/TT = k 
and {z2 — Z3)A{tQ)/7i = k'. The third condition is redundant since Pi + P2 + 
P3 = 1. One may show after some algebra that these conditions are met 
when. 



A{to)/7i = ±no 
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\ 18(n-2)+4(^ 



a = V^2{t)/V2^{t) = -^{n-3) , 

P = Vi3{t)/V23{t) = l , (20) 

where hq is any odd integer. One may then show, 

^[3 + cos(0) +4cos(^/2)] , 

o 

i[3 + cos(^) -4cos(0/2)] , 
8 

■sm'{e/2) , (21) 



Pi = 


|ai(t)p 


P2 = 


|a2(t)P 


P3 = 


\a3{t)\' 



2{n - 2) 

where 6 = e{t) = 2TT[A{t) / A{to)] with A{t) = j/^V{t')dt'. Complete popu- 
lation transfer occurs from state 1 to state 2 at t = to when 6{t) = 27r. As 
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n — > oo these solutions reduce to those of a 2-state atom with Vu = V22, 
where Pi = cos^[A{t)] and P2 = sin^[A(t)]. 

3 Results 

3.1 2-state atom 

In this section we compare full and analytic results for the 2s — 2p transition 
in hydrogen in the two state approximation 

3.1.1 Population leakage 

To obtain our analytic solution to Eq.®, we had to neglect the U2ia2 term 
in the equation 10,2 = uj2i0'2 — x cos(ti;t)ai. This approximation requires that 
^21 X) or (since x/u; = 7r/2), UJ21 ^ u. When UJ21 is finite, the population 
transfer is not complete. The effect of this population leakage can be cal- 
culated by expanding the transition amplitude 02 in a power series in time 
including terms 0021- It is easily shown that the first 0021 term corresponds 
to the second derivative, 02 (t). The difference between the exact and the 
analytic solutions for the transition amplitudes is Aa2(t) ~ \ui2ixt^i and 

AP(t) = |Aa2(t)P^^a;2\xV. (22) 

By the time the occupation probability P2 reaches its first maximum at 
to = T /A = 71 /2u!, the difference becomes 

AP(to) - l&i'^r ■ (23) 

3.1.2 Calculations of 2s — 2p transitions in hydrogen 

Here we compare the full (numerical) calculations and the analytical solution 
for a 2s — 2p transition in hydrogen. The full calculations were done by 
numerically integrating Eqs. (jHl) using a standard fourth order Runge-Kutta 
algorithm. In hydrogen the energy separation of 2s and 2p states (Lamb 
shift) is 4.37 * 10~^eV, while the next available state (3p) is ~ 1.89e\^ away. 
Therefore, one can choose the frequency of the external field uj such that 
LJ21 uj <^ uj2s3p, i-e. both the degenerate state approximation (limit 0021 —>■ 
0) and the two-state model can be used. 
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0.2 0.22 0.24 0.26 0.28 0.3 
Time (in units of T=27:/m) 



Figure 3: Occupation probability P2{t) in the vicinity of to = T/4: analytical 
solution, Eq.Q, and numerical calculations for different values of the ratio 
uj/uj2i (1, 10 and 100). Deviation from the analytic solution decreases as 
iuj2i/ujf (cf. Eqn. (USD) 



As shown in figure El the difference between the full and analytic solutions 
for Pi{t) = 1 — P2{t) in the vicinity oi to = tt/^uj is large when u = 0021- But 
it decreases for u > 0021- The difference is less than 1% for uj/uj2i = 10, and 
0.01% for uj/u2i = 100. Since that difference is proportional to (ci;2i/co')^, one 
can go up in frequency as high as fractions of eV (that brings the difference 
between the approximate and exact solution down to 10^^*^ — 10^^^) and 
yet be far sensibly from the nearest available (2s — 3p) resonant transition 
frequency. Therefore, a radiation field with the wavelength from a few fim 
(and the intensity of the order / cm'^) up to a few cm (and the intensity 

of the order lO'^W/ cm?) can be used for 2s — 2p transitions in hydrogen. The 
probability of the multiphoton excitation to the 3p state (or any other state, 
including continuum) also seems to be small for the range of the external 
field frequencies under consideration. 

Our analytic approximation appears to be valid for 2s — 2p transitions in 
hydrogen over a broad range of field frequencies. 
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3.1.3 Changing the shape of V{t) 

We have noted previously that the choice of the field frequency is a trade-off 
between two competing factors: reducing the population leakage and increas- 
ing the duration of the populated state. If wishes to obtain a long lasting 
populated state with very small leakage, it may be that both requirements 
for the system cannot be met simultaneously for a single frequency radiation 
field. In this case it may be possible to use another form of the interaction 
potential. In particular, one may change the shape of the "fiat" part of the 
probability (cf. figure by using different shapes for the external potential 
Vext{t) ■ For an arbitrary external potential the formulas for the occupation 
probabilities can be written in the form: 

Pi{t) = cos^i f' V2i{r)dT] 
Jo 

P2{t) = sm'[f'v2i{T)dT] . (24) 
Jo 

A Taylor series expansion can be used for choosing the shape of the ex- 
ternal potential. For example, one can use the potential with the first non- 
vanishing derivative of the order higher than four to make the shape of the 
populated state even fiatter. The derivatives can be calculated using the 
general formula for the nth derivative of a composite function JH] , which in 
this case takes the form 

^"-P^-V d-F y' y" y_^ 

dt^ ^i\j\...k\ dy^'hl' 4!^ ■■'M! ^ ' ^ ^ 

where P2{t) = F{y) = siii^ y , y = y{t) = /q V"2i(r) dr , J2 indicates 
summation over all solutions in non-negative integers of the equation i + 
2j + . . . + Ik = n and m = i+ j + ... + k. 

In principle, Eqn. pHjl tells us how to shape V{t) to control the population 
transfer. As one can see from Eq. ()25j) . one may eliminate all terms up to 
the order r'^ by choosing the potential for which all derivatives up to the 
order {k — 1) are zero at the point t = to. That will "fiatten" the shape of 
the populated (or depopulated) state, i.e. allow one to use smaller frequency 
(and, as a result, an increased duration of the state) to achieve the same state 
of population leakage. Therefore, the shape of an external potential can be 
used (along with the choice of the field frequency) for quantum control. For 
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example, in a truly two-state system, choosing the potential of the form 



V(t) = -6{t - to) 



leads to 



Pi(t) = i-e(t-to) 
P2(t) = e(t-to) • 



(26) 



This represents complete and immediate population inversion at t = to- How- 
ever, Fourier transformation of V{t) now contains all frequencies, so one can 
no longer use a two-state approximation for the 2s — 2p transition in hydro- 
gen, since higher states become necessarily involved. 

3.2 3-state atom 

Some allowed values of the action integral, A{to), and the relative interaction 
strength, a, are given in table I. This leads to certain allowed values of 
the action integral A(to) and the relative interaction strength, a (see [12]), 
namely. 



This table includes values of {?T,i,n2}, {no,n'^} and all three {k, k'} cases 
defined in ^2] ■ From more complete numerical output we confirm that rii and 
77-2 each acquire all possible odd integer values, although values of the product 
ni-n2 are restricted. This is consistent with the condition that ni = 2no + n'^, 
and n2 = rio + 2n'^, where Uo and n'^ are arbitrary odd integers. It is also 
evident that the even integer, Ue = Uo + n'^, takes on all even values. One 
may also show algebraically that for each value of the even integer k in case 
i) there are two odd values of an odd integer k in case ii), and vice versa. 
The sets of integers {ni,?T,2}, {no,n'g} and {k,k'} are redundant. The three 
sets of {k,k'} correspond to a single set {r;,i,?T,2}, while the {no,n'^} are in 
one to one correspondence with the {^1,722}. 

Numerical calculations for the time dependence of the populations in 
a degenerate 3-state atom perturbed by external interactions with V{t) = 



(3 = Vn{t)/V2z{t) = ±l . 



(27) 
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ni ■ n2 


rii 


n2 


Tie 


< 




< 


rio 


He 




a 








k 




k 


k' 


-k 


k' 






5 


±1 


±5 


±2 


±3 


±1 


±3 


±1 


±2 


±1.656 


=F2.530 


5 


±5 


±1 


±2 


±1 


±3 


±1 


±3 


±2 


±1.656 


±2.530 


9 


±3 


±3 


±2 


±1 


±1 


±1 


±1 


±2 


±2.221 


0.000 


11 


±1 


±11 


±4 


±7 


±3 


±7 


±3 


±4 


±2.456 


^4.264 


11 


±11 


±1 


±4 


±3 


±7 


±3 


±7 


±4 


±2.456 


±4.264 


17 


±1 


±17 


±6 


±11 


±5 


±11 


±5 


±6 


±3.053 


^5.488 


17 


±17 


±1 


±6 


±5 


±11 


±5 


±11 


±6 


±3.053 


±5.488 


23 


±1 


±23 


±8 


±15 


±7 


±15 


±7 


±8 


±3.551 


^6.487 


23 


±23 


±1 


±8 


±7 


±15 


±7 


±15 


±8 


±3.551 


±6.487 


27 


±3 


±9 


±4 


±5 


±1 


±5 


±1 


±4 


±3.848 


^1.633 


27 


±9 


±3 


±4 


±1 


±5 


±1 


±5 


±4 


±3.848 


±1.633 


29 


±1 


±29 


±10 


±19 


±9 


±19 


±9 


±10 


±3.988 


^7.353 


29 


±29 


±1 


±10 


±9 


±19 


±9 


±19 


±10 


±3.988 


±7.353 


35 


±1 


±35 


±12 


±23 


±11 


±23 


±11 


±12 


±4.381 


^8.128 


35 


±5 


±7 


±4 


±3 


±1 


±3 


±1 


±4 


±4.381 


^0.478 


35 


±7 


±5 


±4 


±1 


±3 


±1 


±3 


±4 


±4.381 


±0.478 


35 


±35 


±1 


±12 


±11 


±23 


±11 


±23 


±12 


±4.381 


±8.128 



Table 1: Some values for the action integral, A{to), and the relative inter- 
action strength, a, allowed for total population transfer. These values are 
found using Eq (j77|) subject to the conditions listed. Here = Uo + n'^, 
Ho = I (2^1 - 77-2) and n'^ = |(2r22 - ^1). 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Time (in units of the field period, T=2k/(x)) 



Figure 4: Occupation probabilities as a function of time, Eq. (jl7|) . The solid 
line denotes P2{t)', the long dash line denotes -Pi(t); and the short dash line 
denotes Psit). In this figure we use the allowed values, a = V12/V23 = 
and A(to) = Jo" V{t')dt' = 2.221, corresponding to rii = 3 and n2 = 3 
{rio = n'^ = 1) in Eq (j7r|) . Complete transfer to state 2 from state 1 occurs at 
t = to = T/4 and again at odd multiples of to- 

Vbcos(ci;t) are presented in figures miHl These results were obtained by using 
a standard fourth order Runge-Kutta numerical integration of Eq()12|) with 
ej = and {3 = 1 for various values of a and Aito). Again complete transfer 
to an initially unoccupied state never occured. However, there were rapid 
oscillations in the populations of all states except near t = T/4 and t = 3T/4, 
where none of the populations oscillated rapidly. This was similar to figure IHl 
This appears to correspond to the onset of complete population transfer, 
which occurs at t = T/4 and odd multiples of T/4, as seen in figures El and El 
Calculations using a few values of a that permit complete transfer to 
state 2 are shown in figures |3 and El Our numerical codes give the same 
results as the analytic expressions of Eq (fT7|) . For n-state systems Eq (j7I|l 
gives population transfers from state 1 to state n. This is shown in figure |31 
where the short dash line corresponds to sum of P3, P4, ■ ■ ■, P„. This provides 
a check that our algebra is correct. We note that < 1 and that the 

^ (ni+n2) — 2 

maximum value of Psit) occurs for rii = n2, where P3 max = |, consistent 
with figure m This corresponds to a = so that direct transitions from state 
1 to state 2 are forbidden. Transfer to state 2 occurs via the intermediate 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Time (in units of the field period, T=2k/(x)) 

Figure 5: Occupation probabilities as a function of time. The solid line 
denotes P2{t)] the long dash line denotes -Pi(t); and short dash line denotes 
Psi^t). In this figure a = —2.530 and A(to) = 1.656, corresponding to ni = 1 
and ^2 = 5 {rio = — 1, n'^ = 3) in Eq(j2Zj). 

state 3. Transfer from state 1 to state 2 and back is complete, and occurs 
periodically. In general a = corresponds to rii = n2 = Suodd, where 
Hodd is any odd integer. This appears to give the simplest condition that 
allows complete population transfer. In this case the action area is A(to) = 
noddle/ V2. 

Calculations for two other values of a that allow complete transfer to state 
2 are shown in figures 0] and El We see that complete transfer occurs twice in 
one period of the oscillating field but that the f reque ncy of the " side bands" 
increases as aA increases. We note that a = i^^^ "^2) becomes either 
large (ni ^ n2, or vice versa), or small (ni ~ ^2) as nin2 increases, while 
A(to) increases as fJnin2/2. When complete population transfer occurs, the 
population lingers in state 2, as seen in figures E] and El It can be shown 
near t = to that 1 — P2{t) varies as [ci;(t — to)]'^• 
Additional control may be achieved by changing the shape of V{t). 
This can be used to control how long the population remains near unity in 
state 2, for example. An interesting example is the case of an ideal, sudden 
'kick' produced by V{t) = AQ6(t — to). If Aq Uodd'^/V^ for example, then 
the mere presence of state 3 allows transfer from state 1 to state 2 without 
any direct transfer from state 1 to state 2. With this ideal 'kick', state 2 is 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Time (in units of the field period, T=2k/(x)) 



Figure 6: Occupation probabilities as a function of time. The solid line shows 
P2{t)', the long dash line shows Pi{t); and the short dash line shows Psit). 
Here a = 8.128 and Aito) = 4.381, corresponding to rii = 35 and n2 = 1 
= 23, n'g = — 11) in Eq(j2Zj). Since these values are allowed, the electron 
population is completely transferred to state 2 from state 1. 

unoccupied before t = to, and fully occupied after t = to. 

To obtain analytic solutions to Eq.Q for a 3-state atom, we imposed 
degeneracy on the 3-state manifold, e.g. by taking Ej = E = 0. \i E ^ 0, 
then E can be removed from Eq(^ by an overall phase transformation, e*^*. 
When the states are not degenerate Uij = {Ei—Ej)/h ^ 0, and the population 
transfer is incomplete. Numerical calculations indicate that the population 
leakage varies as {uJi^jiSf'. The choice of the field frequency, tu, involves a 
trade-off between the duration of time the population remains in state 2 
and the population leakage. The higher the frequency, cj, the smaller the 
population leakage, AP25 but the shorter the duration time, T^, that the 
population remains in state 2. This effect of population leakage is similar 
as that for a 2-state atom PH]- Numerical calculations indicate that this 
population leakage grows rapidly in time. 

4 Discussion 

In this work we have used the 2s — 2j9 transition in hydrogen as an example 
of population transfer. The same approach can be used for any other atomic 
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or molecular system that has a similar pattern of energy states (two states 
located close one to another and far from the other states). As we have 
shown in Eq. ()23p . for an external field with a single frequency tu, the differ- 
ence between the exact calculations and the analytic approximation varies 
as {uj2i/ujY for uj2i/uj <^ 1, and does not depend on the internal structure of 
the atom or molecule. This feature opens the possibility for using different 
systems with different values of the transition frequency, UJ21 and different 
ranges for the field frequency, uo. 

A downside of analyzing more complex atomic or molecular systems is 
that analytic expressions for the orbital functions are not always available. 
Then genetic algorithms (GA) may be used to choose a shape for the ex- 
ternal potential, V{t). The application of GA with active feedback to the 
selective breaking and making chemical bonds in polyatomic molecules has 
been discussed in detail by Rabitz et al ^iTTHIH], and tested experimen- 
tally JHl 1201 )■ ^'^y optimization scheme, the effectiveness of the GA 
increases significantly when the initial value of the parameter close to the 
optimal one is chosen. One might combine our method with GA using ap- 
proximate analytic orbital functions to choose the "starting" form of the 
potential, estimate T, and Per, and then employ the GA scheme. 

For 2s — 2p transitions in hydrogen, sources of the microwave radiation 
with corresponding intensities may now be available, which can be used to 
test our model. Metastable excited H{2s) atoms could interact with a low fre- 
quency (e.g. microwave) radiation. Since the lifetime of the 2p state is much 
shorter than the 2s state, the population of the 2p state may be monitored 
by observing photons emitted in 2p — Is transitions. In these experiments 
the smallness of the Lamb shift requires the use of the temperatures below 
the 1 mK to exclude the thermal transitions between the states. 

In this work we have also considered transitions amplitudes and proba- 
bilities in a 3-state atom with degenerate states. We also assumed that the 
interaction matrix elements, Vij{t) all have a common time dependence. In 
realistic atomic systems, energy states are seldom, if ever, exactly degener- 
ate. Also states outside 3-state manifold usually exist. Consequently, use 
of our results is restricted to external fields with frequencies in the range, 
oJmin < oj < Umax- Here humin is the energy splitting of the nearly degener- 
ate states, and humax is the energy difference between the 3-state manifold 
and the closest state in energy outside the manifold. 

It is instructive to compare our 3-state results to simpler 2-state re- 
sults |141 IT^ . The equations for a 2-state atom may be recovered from 
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our 3-state equations by taking V12 and Vi^ to zero, which corresponds to 
a ^ 00 (with any \ej/a\ < 00), as may be understood from Eq (fT^ . For 
the 2-state atom it has been shown [TH [T3] that aj{t) may be expressed 
in terms of eigenstates Ci{t) via the relation, Ci{t) = Aljjaj(t), where 

A^jj = ^ I ) This matrix may be inverted to give the population 

amplitudes, ai{ t) = ^[y_c i{t) - y+C2{t)] and a2{t) = ^[-Ci{t) + C2(t)], 

where A = 2^Jl + { ^^^'^^ )^ is the determinant of Aiij. A simple solution 

occurs when ei = €2 = e, namely, ai(t) = e~"^*^*) |[e~*'^*^*) + e+*^*^*)] = 
e-^^Mt) cosiA{t)) and a2(t) = e'^^^W Ifg-^^W -e+^^W] = -ie~'^'Mt) sm{A{t)). 
Note that Pi = I — P2. One may then simply determine the conditions under 
which the population of state 2, P2{to), takes on any desired value at time 
to. If ei 7^ €2, it is easily shown that P2(io) < 1/(1 + (^2 - ei)^/4) < 1. 
Thus the diagonal matrix elements of Vij{t) prevent complete transfer to 
the initially unoccupied state when ei and €2 are unequal, and, when they 
are equal, they simply contribute an overall phase. When ei = €2 any 
value of ^2(^0) between and 1 can be found. In particular, P2(^o) = 1 
if A(to) = /o° Vi2{t')dt' = Uodd'^/'i. This allowed value of A(to) for complete 
transfer to state 2 differs from the allowed values given in Eq fl27|l for de- 
generate 3-state atoms. By comparison to a 2-state atom the conditions for 
complete transfer are generally more complex in a 3-state atom, as discussed 
in For example, complete population transfer can occur in a 3-state 

atom when €2 7^ 63. In both cases, however, the duration of time spent in 
the transferred state can be controlled by adjusting the time dependence of 
the external field. Also if the states are not quite degenerate then popula- 
tion leakage occurs. For a harmonic V{t) the population leakage in a 2-state 
atom varies as AP2(^o) — \{jY{^i2/^Y- For 2s — 2p transitions in atomic 
hydrogen complete population control can be nearly achieved in this manner 
using a radiation field with wavelengths (and intensities) ranging from a few 
/xm (with about 10^^ W/cm^) up to a few cm (with about 10^ W/cm^). At 
wavelengths below a few fim coupling to nearby n = 3 atomic states can be 
significant, and this 2-state model breaks down. 

While we have not provided calculations in this work for specific exper- 
iments, some general guidelines for experimental tests and applications are 
evident. First, the model must be valid, so that the frequency of the exter- 
nal interaction is limited by Umin < uj < uomax- If the external interaction 
varies harmonically, V{t) = Vq cos{u!t), then the first line of Eq ()27|) imposes a 
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constraint between Vq and u, namely, Vq/uj = iy^^^Jj^^. For a harmonic in- 
teraction the duration of the time, T^, that the population remains in state 2 
varies inversely with uj; the higher the frequency, u, the smaller the time the 
population remains in state 2. Our degenerate 3-state model can be applied 
to systems with dipole selection rules. A dipole selection rule can correspond 
to a = 0. Transfer from state 1 to state 2 (which is the dipole forbidden 
transition) is complete at time to if the two allowed transition matrix ele- 
ments, V23 and Vi3, are equal in absolute magnitude and A(to) = riodd'^l'^- In- 
some atomic systems state 2 could decay, with a lifetime, T^, to another state 
outside the degenerate manifold and be lost. Such loss can be controlled by 
adjusting uj. Our model also rests on degeneracy. If the three states are 
not degenerate the transfer of population to state 2 is incomplete. We call 
this population leakage. As discussed above this population leakage may be 
minimized by using high frequency external interactions, i.e. uj 3> ujmini but 
at a cost to duration of time, T^, the population remains in state 2. 

Using nearly degenerate states, the duration time, T^, can be further 
controlled J3] by adjusting the shape of V{t). As mentioned in section III 
for example, the transfer is complete, instantaneous and permanent when 
V{t) = riodd-^Sit — to)- Such a quick, hard pulse is called [211 1221 123] a 
'kick'. Two practical limitations on this ideal model are the impossibility of 
producing a signal that varies as 6{t — to), and the existence of an infinitely 
wide spectrum of high frequency components with frequencies uj > ujmax in 
the Fourier spectrum of 6{t — to)- Fortunately these two difficulties can both 
be addressed by using 'kicks' of finite width in time. In some cases it may be 
possible to design a 'kick' so that its duration is short compared to any other 
changes in the system, so that a finite 'kick' may be sensibly represented by 
^{t ~ to)- If) in addition, the energy states outside the (nearly) degenerate 
manifold have a large energy gap hujmax, then it may be possible that uJmin < 
UJ < ujraax and our model may be applicable. Applying Vit) — S(t — to) 
to nearly degenerate atomic systems leads to the 'gedanken' question, what 
happens when one tries to force the transition to occur within a small time 
interval about to in a degenerate quantum system where At is large? 

We note that at sufficiently high uj all bound states in the atomic ele- 
ments become nearly degenerate. Hence, if one can deal with high energy 
continuum states, our approach might be useful in applications involving 
fourth generation synchrotrons that produce intense high frequency fields. 
An extension of this method might also used for control transitions in high 
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Rydberg states pi], including adiabatic rapid passage [^Hllini- KWA is used 
to describe coherent storage of information in photonic states j^Z]. The ap- 
proach developed here may be useful for modeling information transmission 
and storage in atomic states ^| in a new way. 

Features that occur when degenerate states are used and the connection 
to quantum time ordering are discussed elsewhere in this book 

5 Summary 

We have shown that in a ra-state atom with degenerate energies, electron pop- 
ulation is completely transferred via an external interaction at a designated 
time, to, as example we show for n = 2 and n = 3. A new analytical solution 
of two coupled channel equations has been found that enables one to tempo- 
rally control the electron population of 2s and 2p states in hydrogen by using 
a time varying external field. In addition, the population leakage can both be 
easily estimated, and than be further controlled by changing the shape of the 
external field. For 3-state atoms from an initially occupied state (state 1) to 
a designated initially unoccupied state (state 2) under two conditions. The 
first condition for complete transfer is that the ratio of the mat rix el ements of 
the external interaction, l^j(t), satisfy Vi2{t)/V23{t) = a = ±y^^^(ni — 712), 
and Vi3{t) /V23(t) = f3 = ±1. The second condition is that at t = to the 
action area of V{t) satisfy A{to) = V{t')dt' = i-^^f, where we have 
set V23(t) = V{t). Here rii and n2 are integers such that rii = 2no + n'^ 
and n2 = no + 2n'^, where no and n'^ are any odd integers. The duration of 
time the transferred population remains in state 2 can be controlled either by 
varying the frequency, cu, of the external potential, or by varying the shape 
oiVit). 
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